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Abstract 

The stability of a new class of hairy black hole solutions in the coupled system of Einstein- 
Yang-Mills-Higgs is examined, generalising a method suggested by Brodbeck and Strau- 
mann and collaborators, and Volkov and Gal'tsov. The method maps the algebraic system 
of linearised radial perturbations of the various held modes around the black hole solution 
into a coupled system of radial equations of Schrodinger type. No detailed knowledge of 
the black hole solution is required, except from the fact that the boundary conditions at 
the physical space-time boundaries (horizons) must be such so as to guarantee the finite- 
ness of the various expressions involved. In this way, it is demonstrated that the above 
Schrodinger equations have bound states, which implies the instability of the associated 
black hole solution. 
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Coupling gravity to non-linear systems, such as the non-Abelian Yang-Mills the- 
ory, or the non- linear cx-models etc., has led to interesting (classical) solutions with 
particle-like [|TJ or black-hole interpretation 0. The interest in the latter type of so- 
lutions arises mainly from the fact that new types of classical hair have been shown 
to exist, contrary to the no-hair conjecture characterising purely gravitational or 
Abelian black holes ||. This is so, because the no-hair theorems do not involve the 
issue of stability of the solutions in their proof, and therefore in this respect the 
above classical solutions may be considered as explicit counter-examples to these 
theorems. 

In view of this, it is natural to enquire into the stability of the above solutions, 
which would establish their physical significance. It has been shown that most of 
these systems, especially the ones admitting particle-like interpretation, are unstable 
under perturbations of the various field modes 0. For the black hole solutions, a 
corresponding general proof was lacking so far, mainly due to the peculiar behaviour 
of the stability equations on the horizons. In some cases, however, like the Einstein- 
Yang-Mills-Higgs (EYMH) systems with a Higgs triplet, the Einstein-Skyrme (non- 
linear (T-model) system, and the Einstein- Yang-Mills-Dilaton theory (inspired from 
strings), linear stability of the hairy solutions is established |J, although non-linear 
stability remains an unsettled issue. 

An interesting class of classical hairy black holes has been found recently in connec- 
tion with the SU (2)-Einstein-Higgs system, with a Higgs doublet as in the standard 
model M. These black hole solutions resemble the sphaleron solutions in SU(2) 
gauge theory and one would expect them to be unstable for topological reasons. 

Recently, an instability proof of sphaleron solutions for arbitrary gauge groups in 
the EYM system has been given |7], |8j . The method consists of studying linearised 
radial perturbations around an equilibrium solution, whose detailed knowledge is 
not necessary to establish stability. The stability is examined by mapping the sys- 
tem of algebraic equations for the perturbations into a coupled system of differential 
equations of Schrodinger type |7|, || . As in the particle case of ref . , the instabil- 
ity of the solution is established once a bound state in the respective Schrodinger 
equations is found. The latter shows up as an imaginary frequency mode in the 
spectrum, leading to an exponentially growing mode. There is an elegant physical 
interpretation behind this analysis, which is similar to the Cooper pair instability 
of super-conductivity. The gravitational attraction balances the non-Abelian gauge 
field repulsion in the classical solution |]|, but the existence of bound states implies 
imaginary parts in the quantum ground state which lead to instabilities of the so- 
lution, in much the same way as the classical ground state in super-conductivity is 
not the absolute minimum of the free energy. 



However, this method cannot be applied directly to the black hole case, due to 
divergences occuring in some of the expressions involved. This is a result of the 
singular behaviour of the metric function at the physical space-time boundaries 
(horizon) of the black hole. 

It is the purpose of this note to generalise the method of ref. to incorporate 
the black hole solution of the EYMH system of ref. || . By constructing appropriate 
trial linear radial perturbations, following ref. || |J, we show the existence of bound 
states in the spectrum of the coupled Schrodinger equations, and thus the instability 
of the black hole. Detailed knowledge of the black hole solution is not actually 
required, apart from the fact that the existence of an horizon leads to modifications 
of the trial perturbations as compared to those of ref. 0, ||, in order to avoid 
divergences in the respective expressions [f|. 

We start by sketching the basic steps 0, ^| that will lead to a study of the stability 
of a classical solution <p s (x,t) with finite energy in a (generic) classical field theory. 
One considers small perturbations 5<p(x, t) around <p s (x, t), and specifies flU the time- 
dependence as 

5(f)(x,t) = exp{-ittt)ty(x) (1) 

The linearised system (with respect to such perturbations), obtained from the equa- 
tions of motion, can be cast into a Schrodinger eigenvalue problem 

WS> = tt 2 A^ (2) 

where the operators TC, A are assumed independent of the 'frequency' Q. As we 
shall show later on, this is indeed the case of our black hole solution of the EYMH 
system. In that case it will also be shown that 7i is a self-adjoint operator with 
respect to a properly defined inner (scalar) product in the space of functions {^>} 
[0], and the A matrix is positive definite, < >> 0. A criterion for instability 

is the existence of an imaginary frequency mode in (fj) 

n 2 < (3) 

This is usually difficult to solve analytically in realistic models, and usually numer- 
ical calculations are required Q. A less informative method which admits analytic 
treatment has been proposed recently in ref. and we shall follow this for the 
purposes of the present work. The method consists of a variational approach which 
makes use of the following functional defined through (0): 

o2/T . < V\H\V > ... 
n (*) = — TITTF — 4 



with \I/ a trial function. The lowest eigenvalue is known to provide a lower bound 
for this functional. Thus, the criterion of instability, which is equivalent to (|j), in 



this approach reads 

fi 2 (^) < 

< > < oo (5) 

The first of the above conditions implies that the operator 7i is not positive definite, 
and therefore negative eigenvalues do exist. The second condition, on the finiteness 
of the expectation value of the operator A, is required to ensure that \l/ lies in 
the Hilbert space containing the domain of Ti. In certain cases, especially in the 
black hole case, there are divergences due to singular behaviour of modes at, say, 
the horizons, which could spoil these conditions ([5]). The advantage of the above 
variational method lies in the fact that it is an easier task to choose appropriate 
trial functions ^ that satisfy (§) than solving the original eigenvalue problem (|2]). 
In what follows we shall apply this second method to the black hole solution of ref. 

1- 

We start by reviewing the basic formulas for a study of stability issues of spheri- 
cally symmetric black hole solutions of the EYMH system || . The space-time metric 
takes the form |J 

ds 2 = -N(t, r)S 2 {t, r)dt 2 + N~ l dr 2 + r 2 {d6 2 + sin 2 9d(f) 2 ) (6) 

and we assume the following ansatz for the non-abelian gauge potential || |7j 

A = a r r dt + aiT r dr + (u — 1) [r^dO — tq sin 9d(j)] + u[Tgdd + sin 9d(p] (7) 

where u, uj and a^, i = 0, 1 are functions of t, r. The Tj are appropriately normalised 
spherical generators of the SU(2) group in the notation of ref. [0. 

The Higgs doublet assumes the form 

V2 { 0-^3 ; ' V W 

with the Higgs potential 

V($) = j(&$-v 2 ) 2 
where v denotes the v.e.v. of $ in the non-trivial vacuum. 

The quantities uj, (f> satisfy the static field equations 
Nuo" + 
N(j)" + 

where the prime denotes differentiation with respect to r. For later use, we also 
mention that a dot will denote differentiation with respect to t. 
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(10) 



If we choose a gauge in which 5ao = 0, the linearised perturbation equations 
decouple into two sectors . The first consists of the gravitational modes SN, 
8S, 5oj and 5(p and the second of the matter perturbations Sax, Slj and 8ijj. In our 
analysis it will be sufficient to concentrate on the matter perturbations, setting the 
gravitational perturbations 5N and 8S to zero, because an instability will show up 
in this sector of the theory. The equations for the linearised matter perturbations 
take the form |7| 

H^ + A^ = (11) 

with, 




# = | 5u | (12) 

and, 

Nr 2 

A = I 2 | (13) 
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and the components of H are 



H aiai = 2{NS) 2 [u 2 + 
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H n = 2^^ + 2^ 2 (^^±^ + ^A(0 2 -^ 
H ai u> = 2iNS[(p*u) - up*] 

H* ai = 2i[p*NSu + NS(p*u)} (14) 

• 2 
IT 

H ai Tt> = —NS[(p*4>) - 4>p*} 

^2 ^,2 

U^ ai = ip*—NS(f) + i—NS(p*(f)) 
7~(-cjip = Ti-ipQ — —4>NS 2 
where the operator p* is 

p* = -iNs4- (15) 
dr 

Upon specifying the time-dependence (|IJ) 

ty(r,t) = ty(r)e int ; *(r) = I 6u(r) II (16) 




one arrives easily to an eigenvalue problem of the form (0), which can then be 
extended to the variational approach (^j). 



To this end, we choose as trial perturbations the following expressions (c.f. [[7]) 

5ai = u'Z 
6u = (u 2 -l)Z 

81I) = \<j>{u-l)Z 
where Z is a function of r to be determined. 



One may define the inner product 



< ^1$ >= 



1 



(17) 



where is the position of the horizon of the black hole. The operator H. is then 
symmetric with respect to this scalar product. Following ref. |7j], consider the 
expectation value 



, dr NS* 



iVr V) 2 + 2(a; 2 - l) 2 + ^(^ - l) 2 



18) 



which is clearly positive definite for real Z. Its finiteness will be examined later, and 
depends on the choice of the function Z. 



Next, we proceed to the evaluation of the expectation value of the Hamiltonian 
7i QI5D; after a tedious calculation one obtains 

< V\H\V > = / drSZ 2 {-2N(u') 2 + 2P 2 N(u 2 - l) 2 

+\p 2 Nt 2 4> 2 {u - I) 2 - A(u, 2 - I) 2 - U\u - l) 2 } (19) 
+ boundary terms 



where P = ~^?. The boundary terms will be shown to vanish so we omit them in 



the expression (|20). The final result is 
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(20) 



The first of these terms is manifestly negative. To examine the remaining two, we 
introduce the 'tortoise' co-ordinate r* defined by [§] 



dr* 1 
~dr~ ~ Ws 

and define a sequence of functions Z k {r*) by M 



(21) 



r 



Z k {r*) = Z[-) ■ k = l,2,... (22) 



where 



Z(r*) = Z(-r*), 

Z(r*) = 1 for r* e [0, a] 

-D < < 0, for r* e[a,a + 1] 

Z(r*) = for r* > a + 1 (23) 

where a, D are arbitrary positive constants. Then, for each value of k the vacuum 
expectation values of H. and A are finite, < \1 / |7Y| V I / >< oo, and < \I/|j4.|\I/ >< oo, 
with Z = Zk, and all boundary terms vanish. This justifies a posteriori their being 
dropped in eq. (p0|) . The integrands in the second and third terms of eq. (p0|) are 
uniformly convergent and tend to zero as k —>■ oo. Hence, choosing k sufficiently 
large the dominant contribution in (|20|) comes from the first term which is negative. 



This confirms the existence of bound states in the Schrodinger equation (JTTj) , (g), 
and thereby the instability (|5|) of the associated black hole solution of ref. in the 
coupled EYMH system. 

The above analysis reveals the existence of at least one negative odd-parity eigen- 
mode in the spectrum of the EYMH black hole, which implies its instability. The 
exact number of such negative modes is an interesting question and we plan to 
investigate it in the near future. Recently, a method for determining the number 
of the sphaleron-like unstable modes has been applied by Volkov et al. |T(| to the 



gravitating sphaleron case, and one might be able to extend it to the present EYMH 



black hole. According to the analysis of ref. jnj, for EYM black holes, there are 
n - unstable sphaleron-like modes under radial perturbations, where n is the num- 
ber of nodes of the equilibrium solution [|TJ. This number does not depend on the 
details of the equilibrium solution, such as the horizon geometry, size etc. This 
is due to the topological nature of the instabilities. In this respect, we mention 



that an interesting connection could be made with the global analysis of ref. fll 
where catastrophe theory was invoked to provide a way of evaluating the number 
of unstable modes of certain (non-sphaleron) black hole solutions. From the global 



analysis of ref. [ITT] there are other non-sphaleronic types of non-Abelian black holes, 



whose 'high entropy' phase is stable. In our analysis, this would imply an extension 



of the variational approach to incorporate finite temperature effects for the matter 
perturbations in non-sphaleron black holes^. The finite temperature would be a 
result of the existence of the horizon entropy associated with the black hole in a 
semi-classical analysis. It might well be that the number of unstable modes of these 
(non-sphaleron) black holes is somehow affected by the temperature, in the sense 
that above a 'critical' temperature (corresponding to a certain horizon size) the 
bound states of the Schrodinger equation (§) disappear, or their number is reduced. 
This would correspond to the high-entropy 'stable' black holes of ref. [jnj, in the 



sense of the catastrophe theory. At present, such issues remain open. We hope to 
come back to these in the near future. 
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